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Abstract. Double neutron stars are one of the most promizing sources for terrestrial 
gravitational wave interferometers. For actual interferometers and their planned upgrades, the 
probability of having a signal present in the data is small, but as the sensitivity improves, 
the detection rate increases and the waveforms may start to overlap, creating a confusion 
background, ultimately limiting the capabilities of future detectors. The third generation 
Einstein Telescope, with an horizon of z > 1 and very low frequency "seismic wall" may be 
affected by such confusion noise. At a minimum, careful data analysis will be require to separate 
signals which will appear confused. This result should be borne in mind when designing highly 
advanced future instruments. 



1. Introduction 

The coalescence of two neutron stars (BNS), two black holes (BBH) or a neutron star and a 
black hole (NS-BH), are among the most promising sources for ground-based gravitational wave 
detectors due to the huge amount of energy emitted in the last phase of their inspiral trajectory. 
The waveform is very well modeled until the last stable orbit and a detection would provide 
strong constraints on the source parameters. With the third generation interferometer Einstein 
Telescope, the horizon of compact binaries is expected to reach cosmological distances where 
it may become possible to study cosmology. In particular, double neutron stars or neutron 
star-black holes, if associated with an electromagnetic counterpart, may be used as standard 
sirens to constrains dark energy [1]. At such distances however, it is likely that the sources 
create a confusion foreground, where the detection of a single coalescence becomes difficult and 
requires advanced data analysis methods [2]. In the first section, we derive the coalescence rate 
of BNS; in the second section we discuss the different detection regimes and the consequences 
for Einstein Telescope; in the third section we present our conclusions and the work in progress. 

2. Coalescence Rate 

The final merger of a compact binary occurs after two massive stars in a binary system have 
collapsed to form neutron stars or black holes and have inspiralled through the emission of 



gravitational waves. The cosmic coalescence rate is given by: 



oc / ^^P{U)dta with /^°(0) = po (1) 

J l + Zf 

In this expression, pc is the star formation rate in MQMpc~^yr~-^, Zf the redshift at the time of 
formation of the binary system and P{td) the probability distribution of the delay between the 
formation and the coalescence. The coalescence rate per interval of redshift dR° /dz{z) is obtained 
by multiplying by the element of comoving volume dV/dz{z). Here, we assume a distribution 
of the form P{td) oc l/i^ with > tmin and take tmin = 20 Myr as representative of BNSs. 
We use the star formation rates of [31 SI |5l [6] , and consider local rates in the range 0.01 — 10 
Mpc~^Myr^^, with reference models of 1 and 0.4, corresponding to the most current estimates 
derived from statistical studies [T and population synthesis [9l[8]. 

3. Detection Regime 

The contribution of BNS to the instrumental data falls into three statistically very different 
regimes, depending on the duty cycle (or the average number of sources present at the detector 
at the same time): 

A{z)= rr°{z)^{z)dz (2) 
JO dz 

where ^^{z) = P°{z)^ is the coalescence rate per interval of redshift, and t"{z) is the typical 
duration of the inspiral in the detector frequency band, which depends strongly on the low 
frequency limit of the instrument and can last from a few minutes for advanced detectors with 
/l = 10 Hz to a few days for the Einstein Telescope with planned low frequency bound between 
1 - 5 Hz (see Table [l] ) 

(i) Shot noise (A « 1): This case describes when the number of sources is small enough that 
the interval between events is long compared to an individual event's duration. Measured 
waves are separated by long stretches of silence and can be resolved individually. This case 
pertains to instruments that are only sensitive to events at low redshift. 

(ii) Popcorn noise (A ^ 1): As the reach of instruments increases, the time interval between 
events may come closer to the duration of single bursts. Events may sometimes overlap, 
making it difficult to distinguish between them. 

(iii) Gaussian (A >> 1): For instruments with very large reach and excellent low frequency 
sensitivity, the interval between events can be small compared to the duration of an event. 
The signals overlap to create a confusion noise of unresolved sources. 

The limit of the popcorn and the Gaussian regimes, defined by A(z=k) = 1 and A{z^:^:) = 10 
are given in Table [T] [2 for different values of the low frequency bound and local coalescence 
rate. With actual and advanced interferometers, whose horizon is a tens or a hundreds of Mpc, 
we are able to probe only the low duty cycle regime, where sources dont overlap. With the third 
generation Einstein Telescope, on the other hand, the horizon is expected to extend to redshift 
larger than z > 1, and the signal may fall deep inside the confusion regime, especially between 
1 — 5 Hz, where it can last for a few days (see Table [T]). 

We find that if the low frequency bound is at /l = 1 Hz, the transition to a Gaussian 
background at z^:^: is well within the ET horizon; 'd = 5 Hz, it most likely occurs at 
z^:^: ~ 0.6 — 1.2, but can fall beyond the detection horizon if p° < 0.15 Myr~^ Mpc~^. The 
transition to a popcorn background at z^,, on the other hand, always occur before the detection 
horizon, and more likely around ~ 0.25 — 0.4, unless for the most pessimistic coalescence rates 
(p° < 0.015 Myr^i Mpc'^). 
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Fi gure 1. Evolution with the detector horizon Zmax of a GW time serie from a simulated 
population of BNSs. Here we have assumed detector a low frequency bound = 10 Hz. 
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Figure 2. Evolution with the detector low frequency bound, of a GW time serie of the from a 
simulated population of BNSs. Here we have assumed a detector horizon Zmax = 0.5. 



Table 1. Threshold between resolved and unresolved BNSs for different estimates of the source 
rate /9° and detector lower frequency bound fi- No value means that the number of sources at 
the detector is always < 1 or < 10. tq refers to the typical duration of a source located at z = 0. 
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4. Conclusions and future work 

With the first and second generations of gravitational-wave detectors the chance of detection 
is limited by the instrumental noise, but as the sensitivity and the number of sources increase, 
the waveforms may start to overlap, creating a confusion background, ultimately limiting the 
capabilities of these instruments. Experience from the Mock LISA Data Challenges [lOj and 
ideas developed for the Big Bang Observatory prove that disentangling multiple signals in 
a gravitational-wave detector's data stream is certainly possible. We have developed a Monte 
Carlo simulation code to generate the signal from populations of extra-galactic binaries in the 
output of detectors, which can be used to study the statistical properties of the background 
and investigate advanced data analysis strategies. The evolution of the signal with the detector 
horizon and with the low frequency bound is illustrated by the simulated time series in Fig. [T] 
and Fig. |2j One can imagine applying a high pass filter to identify and subtract the sources, 
iteratively from lower to higher redshifts, from the integrated signal. This work is in progress 
and will be reported in a future paper. 
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